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Abstract 

Nonlocality is a property of paramount importance both conceptually and computationally ex- 
hibited by quantum systems, which has no classical counterpart. Conceptually, it is important 
because it implies that the evolving system has information on what happens at any space point 
and time. Computationally, because such a knowledge makes any calculation intractable as the 
number of degrees of freedom involved increases beyond a few of them. Bohmian mechanics, with 
its trajectory-based formalism in real configuration space, can help to better understand nonlo- 
cality. A detailed analysis of how nonlocal information is transmitted to quantum trajectories 
in simple systems (free particle and harmonic oscillator) turns out to be very interesting when 
compared to analogous systems in classical mechanics. 
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I. INTRODUCTION 



Since the early steps of quantum theory there has always been much interest in estab- 
lishing a connection between this theory and classical mechanics. A lot of work has been 
devoted to derive quantization rules based on the correspondence principle 111, to propose 

□ n n n 

and improve semiclassical approximations |2J, y, m, |5j , or to develop quantum-classical stud 
ies by using the Liouvillian formulation of quantum mechanics [6 



13, 

approximations where quantum-classical mixed dynamics are involved 



or some sort of 
10j. However, there 



is an alternative pathway consisting of exploiting the properties offered by the so-called 
quantum Hamilton- Jacobi equation, as done by Leacock and Padgett [ll| in 1983. This 
equation can be derived from the time-dependent Schrodinger equation by means of a sim- 
ple transformation which allows to express the wave function in terms of a complex phase 
(which also contains all the information that carries the wave function). The formulation 
based on the quantum Hamilton- Jacobi equation was formerly aimed to obtain bound-state 
energy levels of quantum systems with no need to solve the equation of motion involved in 
the calculation of the system wave function. 

Leacock and Padgett formulation leads straightforwardly to quantum trajectory ap- 
proaches, which have also been an alternative to understand the relationship between quan- 
tum and classical mechanics. Bohmian mechanics, which was proposed before, is also one 
of these approaches. Nonetheless, it presents an inconvenient when dealing with (energy) 
eigenstates. These states are associated with zero velocity fields, and therefore particles 



Faraggi and Matone 
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remain steady when described by them [12|. To solve this problem, Floyd [13j, ll4( and 



16| developed quantum Hamilton- Jacobi like formulations start- 



ing from (real) bipolar ansatzs, though they did not claim full equivalence with standard 
quantum mechanics regarding their predictions. This problem was also treated by John 



17| | by means of a theory based on complex quantum trajectories derived from the ansatz 



used by Leacock and Padgett, and that was applied to some simple analytical cases, such 

ator and the step barrier. This trajectory-based approach has been 



as the harmonic oscil 
used recently 



20 



21] to propose different numerical algorithms that allow to obtain 



transition probabilities and bound-state energies. 

When quantum trajectory-based approaches are used, a full localization underlies the 
particle dynamics since we are dealing with well-defined trajectories in (either real or com- 
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plex) configuration space. Then, one might have the wrong impression that nonlocal effects 
are somewhat washed out. Remember that nonlocality is present in so remarkable phe- 
nomena such as entanglement, which is the cornerstone of quantum information theory and 
quantum computation, for instance. Hence nonlocality becomes an important issue when 
working with quantum trajectories (mainly when proposing alternative methods to solve 
quantum problems by means of them), because not only it does not disappear at all, but 
plays a key role in the particular topology displayed by such trajectories. Independently 
of connotations coming from entanglement and quantum information theory, nonlocality is 
present even when considering the time evolution of a single particle. This can be nicely 
seen by means of Bohmian mechanics, where quantum trajectories are defined in the (real) 
configuration space where experiments are carried out. In order to reproduce experimental 
outcomes, one needs to collect many particles which (we assume) describe certain trajectories 
within the Bohmian approach. All these particles are independent (i.e., there is no physical 
interaction among them) and reach the detector (whatever it is) one by one. Despite their 
independence, all of them are "guided" by identical initial wave functions (or wave fields); 
the corresponding trajectories are then obtained from the time evolution of the wave field, 
and the outcomes from the asymptotic limits in space and time of p(x,t) = \^(x, t)\ 2 (we 
can consider one dimension in this description without loss of generality). The effects of \1/ 
on the particle motion can be then considered as a manifestation of how nonlocality arises 
in quantum mechanics. 

In this Letter a discussion on nonlocality is carried out in terms of the WKB-like formu- 

nn 

lation of Bohmian mechanics (WKB-BM), proposed by Sanz et al. [51. 122l. 



23l | to understand 



the quantum-classical correspondence. This correspondence has been recently considered 



in 



;he study of different problems appearing in realistic atom-surface scattering problems 
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26|. 



II. WKB FORMULATION OF BOHMIAN MECHANICS 
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The starting point of the WKB-BM approach consists of using the general ansatz 

V(x,t) = e lS ^' h } (1) 
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where S is a complex (phase) function. This simple relation allows to express the time- 
dependent Schrodigner equation as a quantum Hamilton- Jacobi equation, 

dS (V5) 2 th~- 

Equivalently, applying the V operator, this equation can be expressed in its hydro dynamical 
form, 

— = + — V 2 v, (3) 

at m 2m 

where v = VS/m is a complex velocity field and d/dt = d/dt + vV is the associated La- 
grangian operator. This equation is the starting point of different numerical algorithms 



proposed in the literature [20|, |27| to deal with a large variety of quantum problems, partic- 



ularly high- dimensional ones, where standard quantum mechanics is difficult to apply. 

After Eq. (j2J), the WKB-BM approach considers an expansion of S as a power series of 
h/i, 

S(x,t) = Y,(-) S n (x,t), (4) 

n=0 V % J 

where the functions S n are real (as seen below, they satisfy a hierarchy of real coupled 
equations). We are interested in knowing how and when nonlocality appears, and therefore 
whether a motion can or cannot be regarded as classical. Therefore, the functions S n have 
to obey the usual conditions of the standard WKB approximation [28]. Substituting (J4j) 
into Eq. (j2J), the latter can be conveniently decoupled in powers of h/i to yield 

dS , (V£ ) 2 



dt 2m 



+ V = (5) 



for the zeroth order (i.e., the classical Hamilton- Jacobi equation, with So being the classical 
action), and 

n _ 

+ ^7 E V5 fc V5 n _ fc + — V 2 ^x = (6) 



dt ' 2m ^— ' " " 2m 

k=0 



for higher orders, n > 1. Eq. ([6]) together with Eq. (JSJ) provides an iterative scheme to 
compute higher orders in the expansion of S. In particular, when the series is truncated 
at n — 1 one gets the well-known WKB approximation. This scheme is not based on the 
propagation of individual trajectories, but on the propagation of functions displaying the 
global behavior of ensembles of particles. In this sense, the S n functions can be considered 
as fields. 
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Once every S n is known, the series (JH) will also be known. One can then compute the 
real fields R and S that define the Bohmian (polar) ansatz 

V{x,t) = R{x,t)e lS{x ' t)/h (7) 

in terms of the S n functions by identifying Eqs. ([!]) and ((7j) to give 



R = exp 



^(-l) n h 2n S 2n+1 



n=0 



(8a) 



S = ^(-l)«n 2 ™ S 2n . (8b) 

n=0 

We would like to stress that although S is given in terms of even powers of h, eventually 
one could observe constant or time-dependent terms depending linearly on H (see below). 
Note also that the wave function is always well defined except an overall phase that does 
not depend on the coordinate space, which is not taken into account in (DO). By substituting 
Eq. ( I8bl) into the Bohmian velocity field, 

x = v = , (9) 

m 



we reach 



71=1 



1 OO 

+ - y^(-i) n h 2n vs 2n , (io) 



where x^ cV) = VSo/m is the classical law of motion. From Eq. (TTUT) one can therefore interpret 
quantum trajectories as classical trajectories "dressed" with a series of terms coming from 
quantum interference. This result makes apparent the important difference between both 
types of trajectories, classical and Bohmian. Within the WKB-BM approach, nonlocal 
effects are contained in the second term of the r.h.s. of Eq. fflOl) and start appearing at very 
small times. 



III. APPLICATIONS OF THE WKB-BM APPROACH 

In the light of the WKB-BM approach the effects of the nonlocality are nicely illustrated 
by means of the following well-known examples. First, we consider the case of the free 
propagation of a Gaussian wave packet, 

/ i \ 1/4 

\£^ x t) = I J e -(x-v t) 2 /45 t aa+ipo(x-v t)/h+iEt/h ^ 

\2ira t J 



where the time complex spreading is given by at = cr (l + 2ftt/2mcx 2 ), the initial velocity 
by vo = Po/m = (p/m), and the initial energy by E = Po/m = (H). The corresponding 
(Bohmian) action is then 

ni \ h ( ftt \ r-, h 2 t 2 

S(x,t) = -- arctan ^ + £/t + + 3 5-0 - f t) , (12) 

with cr t = (To[l + (ftt/2mcr 2 ) 2 ] 1//2 being the time- dependent real spreading. For ft small (or, in 
general, htj1ma\ small), note that S can be expressed, effectively, as a series of even powers 
of h when the arctangent function and the spreading of the wave packet are expanded in a 
Taylor expansion. 

Introducing (fl2l) into Eq. (jUj), we obtain the exact expression for the quantum trajectories, 

x(t) = v t H — -x , (13) 

where the initial condition is x(0) = xq. If we assume htj2ma\ < 1, (TIB"]) can be expressed 
as 

2n 



x(t) = Xq 



V^V ^fn i) (2n-3)!! / ftt \" n 

n=l v 07 

As expected, the two first terms in the r.h.s. of this expression are exactly the same that one 
would expect from a classical trajectory; note that the first term arises precisely when there is 
no spreading. However, the presence of the ft-dependent (third) term makes that a divergence 
with respect to the classical motion be observable: it leads to a hyperbolic spreading of the 
trajectories. This is a very remarkable effect where the nonlocality, as defined above, plays 
an important role; in order to avoid crossing among trajectories, those with the outmost 
initial conditions (with respect to Xq = 0) will spread at a faster rate than the innermost 
ones. This information is contained in the spreading ratio a t /ao, which provides global 
information (i.e., it transmits the information carried by the wave function) on how the full 
velocity field has to be at any time. In other words, in this simple example the second term 
of Eq. ([TBI is responsible for nonlocal effects. Nevertheless, note that for relatively small 
h (or, equivalently, very short timescales, wide wave packets, and/or massive particles), 
the quantum trajectories will basically follow a similar behavior as classical particles, this 
explaining the good agreement between quantum mechanics and semiclassical approaches in 
this case. 

The general trend described above is illustrated in Fig. 1, where a sample of Bohmian 
and classical trajectories with equal initial conditions are shown. Observe that after a very 




FIG. 1: Bohmian trajectories (solid lines) corresponding to a free Gaussian wave packet and 
obtained by means of Eq. (fT3j) at: (a) short timescales (region around tt/2mag < 1) and (b) 
asymptotic times (ht/2maQ ^> 1). To compare, the corresponding classical trajectories have also 
been represented with dashed lines. In part (b), dashed lines indicate the classical- like asymptotes 
described by Eq. (fT5|) . 



short time, Bohmian trajectories start displaying a completely different behavior than their 
classical counterparts. Moreover, note also that at htj2mo\ ~ 1 the series given by Eq. (|T4|) 
breaks down and can no longer be used to describe the long-time behavior of the Bohmian 
trajectories. It is then interesting to consider the asymptotic limit of Eq. (|T3|) . given by 

hxo 



x(t) 



(15) 



This result shows us that asymptotically (within the so-called Fraunhofer regime of motion 



22 



23, 



29 



301 ]) quantum trajectories (as classical trajectories) describe a uniform rectilinear 
motion. However, the corresponding quantum constant velocity given by v' Q = vo + hxo/2m<jQ 
contains a residual term coming from the nonlocality (the velocity of each particle depends 
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on its corresponding initial position). As discussed elsewhere [22j, |23|, |2J, |25|, |26|, [29|, |30| . 
the nonlocal and context-dependence information transmitted to the quantum trajectories 
is through the quantum potential, which is related to the curvature of the wave function 
and therefore is responsible for the residual contribution of the constant velocity. 

Next example is the harmonic oscillator, which is also worth discussing. In this case, the 
time-evolution of a Gaussian wave packet centered around x = a is described by 

/ i \ 1/4 

V]/^ ^ j j g— (x— a cos a;t) 2 /4<TQ —iuit/2— imuj{ixa siiiuit— a? sin 2iot)/ih (16) 

where a 2 = h/2muj. Here, the quantum action is 

| TflUJ 

S(x,t) = — — Tkut — (Axasinut — a 2 sin 2ut) , (17) 

and therefore the quantum trajectory is given by 

x(t) = (xq — a) + a cos cut, (18) 

which satisfy the motion equation x + u 2 x = uj 2 {xq — a). As clearly seen from fjlSp . particles 
will display oscillations parallel to both the classical trajectory (which oscillates around 
x = a with frequency u) and other quantum trajectories [i.e., given two particles, labeled as 
1 and 2, their respective trajectories will satisfy the property x^(t)—x^(t) = x^ — Xq L) \. To 
understand how nonlocality operates here, note that the wave packet does not spread with 
time but remains the same, with its center following the path tracked by the corresponding 
classical trajectory. This implies that the quantum motion is constrained. Thus, since 
trajectories cannot cross, the only possibility for their topology is to be the same as the 
classical one (different initial positions will give different parallel trajectories). That is, the 
nonlocality manifests in the distribution of initial conditions, but not in the particular value 
of h, since out = ht/2maQ can acquire any value. 



IV. CONCLUSIONS 



As shown here, nonlocality is strongly related to the fact that, at any time, the evolution 
of a (quantum) system strongly depends on the full configuration of the real (coordinate) 
space, rather than the features presented by the particular formulation used to describe it. 
In other words, the topology of the trajectories is strongly affected by the information carried 
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in the full wave function (nonlocality could be thus described saying that at every time any 
particle has information about the whole system configuration, unless h = 0). Therefore, 
the nonlocal behavior arises from having a complete information encoded within the wave 
function, which is transmitted to the (quantum) particles and indicates how they should 
evolve according to certain rules as illustrated with the above two examples. Due to the fact 
that the classical limit h~ — > is carried out analytically, one can then observe at each time 
step the classical and quantum behaviors. Of course, as the system becomes more complex, 
the classical limit is not carried out analytically and nonlocal effects are not easily detected. 
When mixed dynamics 10j are proposed to solve problems with high dimensionality, a 
detail analysis of the nonlocality should also be carried out in order to give a complete 
interpretation of the results in terms of trajectories because residual contributions, as those 
seen in Eq. ffl5|) . can be masked by a fully classical analysis. 

On the other hand, we would like to emphasize the insight that Bohmian mechanics 
renders when applied to the study of physical systems. In general, trajectory schemes can 
be of great help to develop efficient computational algorithms to solve quatum problems. 
However, at some point, one should also balance this efficiency with the insight that they 
can provide. Here we have seen that a WKB-BM scheme allows a nice comparison between 
Bohmian and classical trajectories and, at the same time, constitutes a method to tackle 
a quantum problem at different levels of approximation depending on the (relative) value 
of h. In particular, this scheme is useful when trying to compare real Bohmian trajectories 
with trajectories obtained from the WKB wave function. Generally, when using the WKB 
approximation, one only calculates the first terms in the expansion 2J, |25|, |26| . So and 
Si, which already contain nonlocality, as seen above. Nevertheless, sometimes it could be 
necessary to go to higher orders in semiclassical approximations as happens, for example, 
for deep tunneling where one needs more terms in order to account for the longer tunneling 



times involved in the process 



3l\ . Finally, we would like to emphasize that due to the large 



number of works in WKB and WKB-like formalisms, the Bohmian analysis as proposed here 
could be of great help to provide a deeper inside in nonlocal effects. 



9 



Acknowledgements 



Support from the Spanish Ministry of Education and Science under the project No. 
FIS2004-02461 is acknowledged. A.S. Sanz would also like to thank the Spanish Ministry of 
Education and Science for a "Juan de la Cierva" Contract. 



[1] P. Brumer, J. Gong, Phys. Rev. A 73 (2006) 052109. 
[2] F. Grossmann, Comments At. Mol. Phys. 34 (1999) 141. 
[3] W.H. Miller, J. Phys. Chem. A 105 (2001) 2942. 
[4] S. Zhang, E. Pollak, Phys. Rev. Lett. 91 (2003) 190201. 

[5] R. Guantes, A.S. Sanz, J. Margalef-Roig, S. Miret-Artes, Surf. Sci. Rep. 53 (2004) 199. 

[6] C. Jaffe, P. Brumer, J. Phys. Chem. 88 (1984) 4829. 

[7] C. Jaffe, P. Brumer, J. Chem. Phys. 82 (1985) 2330. 

[8] J. Wilkie, P. Brumer, Phys. Rev. A 55 (1997) 27. 

[9] J. Wilkie, P. Brumer, Phys. Rev. A 55 (1997) 43. 
[10] E. Gindersperger, C. Meier, J. A. Beswick, J. Chem. Phys. 113 (2000) 9369. 
[11] R.A. Leacock, M.J. Padgett, Phys. Rev. D 28 (1983) 2491. 

[12] P.R. Holland, The Quantum Theory of Motion, Cambridge University Press, Cambridge, 1993. 

[13] E.R. Floyd, Phys. Rev. D 26 (1982) 1339. 

[14] E.R. Floyd, Int. J. Mod. Phys. A 14 (1999) 1111. 

[15] A.E. Faraggi, M. Matone, Phys. Lett. A 249 (1998) 180. 

[16] A.E. Faraggi, M. Matone, Int. J. Mod. Phys. A 15 (2000) 1869. 

[17] M.V. John, Found. Phys. Lett. 15 (2002) 329. 

[18] C.-C. Chou, R.E. Wyatt, J. Chem. Phys. 125 (2006) 174103. 

[19] C.-C. Chou, R.E. Wyatt, Phys. Rev. E 74 (2006) 066702. 

[20] Y. Goldfarb, I. Degani, D.J. Tannor, J. Chem. Phys. 125 (2006) 231103. 

[21] A.S. Sanz, S. Miret-Artes, J. Chem. Phys., submitted for publication. 

[22] A.S. Sanz, F. Borondo, S. Miret-Artes, J. Phys.: Condens. Matter 14 (2002) 6109. 

[23] A.S. Sanz, Ph.D. Thesis, Universidad Autonoma de Madrid, Madrid, 2003. 

[24] A.S. Sanz, F. Borondo, S. Miret-Artes, Europhys. Lett. 55 (2001) 303. 



10 



[25] A.S. Sanz, F. Borondo, S. Miret-Artes, J. Chem. Phys. 120 (2004) 8794. 

[26] A.S. Sanz, S. Miret-Artes, J. Chem. Phys. 122 (2005) 014702. 

[27] R.E. Wyatt, Quantum Dynamics with Trajectories, Springer, New York, 2005. 

[28] K. Gottfried, Quantum Mechanics, W.A. Benjamin, New York, 1966, Vol. 1. 

[29] A.S. Sanz, F. Borondo, S. Miret-Artes, Phys. Rev. B 61 (2000) 7743. 

[30] A.S. Sanz, S. Miret-Artes, J. Chem. Phys. 126 (2007) 234106. 

[31] D.H. Zang, E. Pollak, Phys. Rev. Lett. 93 (2004) 140401. 



11 



